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Abstract 

We present a simple method to derive the semiclassical equations of motion for 
a spinning particle in a gravitational field. We investigate the cases of classical, 
rotating particles, i.e. the so-called pole-dipole particles, as well as particles with 
an additional intrinsic spin. We show that, starting with a simple Lagrangian, 
one can derive equations for the spin evolution and momentum propagation in the 
framework of metric theories of gravity (general relativity) and in theories based on 
a Riemann-Cartan geometry (Poincare gauge theory), without explicitly referring 
to matter current densities (spin and stress-energy). Our results agree with those 
derived from the multipole expansion of the current densities by the conventional 
Papapetrou method and from the WKB analysis for elementary particles. 

PACS number: 04.50+h, 04.20.Fy 

1 Introduction 



It is of interest to have an alternative method of deriving the Mathisson-Papapetrou 
equations (see pQ and [2]) not only because the conventional procedure is rather lengthy, 
but more importantly, because one might get a better insight into the interpretation of 
the various quantities involved (see also jH] for a detailed analysis of the classical method). 
For instance, in Riemann-Cartan space, there exists a large variation of possibilities of 
how to define the momentum vector which eventually lead to differences in the final 
equations (see [I] and [3], especially the remarks at the end of 0). Furthermore, if we 
consider a specific solution for the metric and the torsion fields, it would be convenient to 
have a Lagrangian into which we can plug the results and derive the equations of motion 
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by varying the Lagrangian, instead of plugging the solution directly into the equations 
of motion. For instance, in general relativity, one usually substitutes the Schwarzschild 
metric into the Lagrangian L = mgikU l u k and not into the geodesic equation, which 
would require to evaluate the complete Christoffel symbols. In this way, one can take full 
advantage of the symmetries of the problem. The aim of this article is therefore to find 
suitable generalizations of L = mgikU l u k that allow for the description of particles with 
spin (classical and intrinsic) moving in a Riemann-Cartan geometry. 

To avoid confusion, let us make a remark on the language we will use throughout this 
article. We will use the expression classical spin of a body to denote its orbital angular 
momentum in its (center of mass) rest frame. The word classical is used to distinguish it 
from the intrinsic spin which is of quantum mechanical origin, and the word spin is used 
to emphasize that we are always talking about the angular momentum in the rest frame. 
When talking about intrinsic spin, we either refer to the intrinsic spin of an elementary 
particle, or to the total intrinsic spin of a macroscopic body as a result of the polarization 
of its constituent particles, for instance in a ferromagnet or in a neutron star. The word 
spin may refer to both classical and intrinsic spin. 

Several attempts have been made for a Lagrangian description in the past and results 
were achieved not only in general relativity for the classical multipole particle [6 , but 
also for point particles with intrinsic spin in a Riemannian space [7], [H|. We will review 
these cases, discuss differences in our derivation method and generalize to particles with 
intrinsic spin and dipole moment (classical spin) in a general Riemann-Cartan space. 

Let us use the example of a test particle in an electromagnetic field to illustrate how 
we will manage to achieve equations for the position as well as for the spin evolution, from 
a single Lagrangian, without referring to the field equations or the Bianchi identities. 

The Lagrangian of a charged particle reads 

L = —g^u* - eAiu\ 

The charge e is of course the integrated current j m and its conservation law e = (the 
dot denotes the time derivative with respect to proper time r) can be derived from the 
relation j m m = which is a consequence of the field equations. However, we do not want 
to refer to the field equations at all and we will ignore the relation e = for the moment. 
The equations of motions are then derived in the form 

mDu 1 = eF\u k +A { —. 

k dr 

The symbol D denotes the covariant derivative with respect to the proper time r. Without 
knowing the Maxwell equations, we certainly know that they are gauge invariant with 
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respect to A rn — > A m + / >m and consequently this should be the case for the equations 
of motion too. We therefore have to require e = 0. In this way, we are let to the charge 
conservation equation as well as to the correct equation of motion. How do we know 
that the field equations will be gauge invariant? Well, the action for the Maxwell field is 
constructed under exactly this requirement. The gauge invariance is more fundamental 
then the field equations themselves. It incidentally leads to the unique Lagrangian F 2 . 
(This will not be the case for the local SO(3, 1) gauge invariance of Poincare gauge theory, 
which allows for a whole family of Lagrangians.) 

In the same way, we will derive the equations for the evolution of the intrinsic spin 
and/or classical spin by investigating the covariance properties of the momentum equation 
with respect to coordinate transformations and/or Lorentz gauge transformations. 

In the next section, we derive in this way the Mathisson-Papapetrou equations for the 
pole-dipole particle in general relativity. In section 3, we generalize the result to Riemann- 
Cartan geometry and finally in sections 4 and 5 we investigate the case of particles with 
intrinsic spin and of macroscopic spin polarized bodies. To get more confidence in our 
procedure, we will re-derive in section 6 the equations of motion for the case of the 
Dirac particle using an alternative method. In section 7 we consider as an illustrative 
example the precession equation of a spinning particle in the torsion field arising in Brans- 
Dicke theory generalized to Riemann-Cartan space. Finally, in section 8, we analyze the 
symmetry properties of the Lagrangian itself, as opposed to the equations of motion. 

2 The classical Mathisson-Papapetrou equations 

In this section, we consider a classical (i.e. without intrinsic spin) pole-dipole particle 
moving in a curved, Riemannian spacetime. 

Instead of postulating a Lagrangian, we will try to derive it from simple considerations. 
The following manipulations, although not mathematically rigorous, will provide us with 
a strong motivation for the Lagrangian (8) which will be our starting point. 

In general relativity, the geodesic equations for a body of mass m are derived from the 
following Lagrangian 

m i k 
L = —9ikU u . 

First, consider the body to be constituted from freely moving test particles, moving all in 
a given exterior field. Then, the Lagrangian can be written as 

L = ^^g ik (x {l) )u\ l) tf l) . (1) 
i 
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However, in order to describe an extended, rotating body, the coordinates of the 
different constituents will have to fulfill certain constraints, due to the rigid structure of 
the body. We now follow the method that is usually used on post-Newtonian Lagrangians 
(see jH] for instance) and also in classical mechanics, but we will apply it directly to the 
covariant form (1). Let X % denote the center of mass of the body. (We avoid the difficulties 
and ambiguities concerning the definition of the center of mass in general relativity. We 
just suppose that there exists a vector X 1 whose spatial components X a , a = 1, 2, 3 agree 
in some limit with the Newtonian center of mass.) 

Then, for the coordinate of a mass element of the body we can write (we omit the 
index denoting the mass element) 

x* = X* + p\ (2) 

and for the velocity of the same mass element 

u i = V i + u\p k , (3) 

with antisymmetric uj tk and V 1 = dX*/dr. Just as in the original work of Papapetrou j2], 
we suppose p° = in a certain reference frame. If we set in the same frame a/ = 0, the 
quantity u %k is clearly the angular velocity in four dimensional form, i.e. we have for the 
spatial part and the three velocity reads u = V + tu x p . Consequently, uj 

is related to the angular momentum (the classical spin) through 

S lk = Iu ik , (4) 

where / is the moment of inertia of the body, defined as in classical mechanics, i.e. we 
have in the spherical case I = I xx = I yy = I zz , with = J (—g^p 2 — p^p u )dm. 
We will use the three dimensional relation 

\l<T = - J P^p'dm. (5) 

This relation can be found in [1^] (§106, problem 4) and also in 9J. It is easily proved 
by taking its trace and using the fact that I xx = I yy = I zz . (In the expression for the 
moment of inertia, the metric can be considered to be independent or p, since the integrant 
is already of order p 2 . See expansion below.) Note that (5) is essentially a Newtonian 
order relation. There is no univoque generalization of the concept of moment of inertia 
in the framework of general relativity. 

Since by definition, X % is the particle's center of mass, we have 

J p ! dm = 0. (6) 
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This equation, again, can only be correctly understood at the Newtonian level. We do 
not specify the integration meter dm in general, but just assume that it has the correct 
Newtonian limit. We now write (1) in the form L = J ^g,ik{x)u l u k dm and expand the 
metric gi k (x) = gi k {X) + g i k i i(X)p l . By using (3)- (6), and retaining only the first order 
metric perturbations and terms at most linear in u %k we find 

777 I 

L = -g ik V l V k - -g ik , m S im V k , (7) 

where V 1 = dX*/dr. Note that the terms linear in p do not contribute because of (6). 
Since p % has been eliminated, we can now leave the frame where p° = u) t0 = and consider 
S vm as a four dimensional antisymmetric tensor. 

Using the symmetry properties of the Christoffel symbols T l km and of S im , we can 
finally write (7) in the simple form (we change the notation from V % to u % ) 

777 1 -"• 

L = -g lk u*u k --T km S k u m . (8) 

This Lagrangian (generalized to higher order multipole terms) has been used in jH], start- 
ing from Newtonian quantities and looking for four dimensional generalizations. A La- 
grangian approach for spinning media can also be found in [TTj . 

Let us emphasize once again that the whole dervivation that led to (8) is essentially 
based on Newtonian physics and can therefore not be used to conclude, or even to prove, 
that (8) is the correct Lagrangian for the description of the dipole particle. The arguments 
make however clear that (8) is probably a good candidate to start with, because by 
construction, it will certainly lead to equations that possess the correct Newtonian limit. 
Beyond that, we can only compare the equations of motion dervived from (8) with those 
derived by other means (as the Papapetrou method) or directly with experiment. In this 
sense, we simply postulate the Lagrangian (8) and forget about its origin. 

Since S ik is an integrated quantity, it is a function of proper time only, without explicit 
x l -dependence. Therefore, the Euler-Lagrange equations are found to be 

mDu m = ~(f l ki , m - f L,JS, V + \? k Ji k - (9) 

Next, the authors of |E] proceed as follows. In a coordinate system where f \ k = 0, 
equation (9) can be written as 

mt)u m = -^R l kmt S l k u\ (10) 
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and since this is a tensor equation, it has to be valid in every coordinate system. In 
this way, starting with a more general Lagrangian, they were able to derive equations of 
motion very effectively, including up to octupole moment contributions. Apart from (10), 
the spin evolution equation was written down as covariant generalization of the Newtonian 
precession equation. In brief, this method consists in looking for covariant equations that 
have the correct Newtonian limit. In a Riemannian theory, this is a very powerful tool 
and is infallible in the sense that there is only one covariant generalization to a given 
Newtonian equation, as long as we suppose that the minimal coupling prescription is 
valid (that is, matter fields should not couple directly to curvature). 

Here, we will follow a different procedure. Since we will later on generalize our ge- 
ometry to Riemann-Cartan space, there will in general be more than one way to write 
a certain Newtonian equation in a covariant form. For instance, both the geodesic and 
the autoparallel equation, although they differ in a general Riemann-Cartan geometry 
and although they are both covariant, reduce to the correct, free particle motion in the 
absence of gravitational fields (curvature and torsion). Therefore, the mere knowledge of 
the Newtonian limit is not enough to fix the general form of an equation. 

We proceed as in the electromagnetic analogy in the introduction. We claim that (9) 
is already covariant as it stands and look at the consequences of this claim. If we complete 
the first two terms of the r.h.s. of (9) to form the curvature tensor, we find 

mDu. 



It is now easy to recognize the covariant derivative f)S t k in the second row 
have 

mbu rn = —Bf^S^u* + ^DS,*. 
For this equation to be covariant, we have to require in addition 

DS ik = 0. (12) 

In this way, we finally get equations (10) and (12) from the Lagrangian (8). They are in 
agreement with the Mathisson-Papapetrou equations in the required order of precision, 
i.e. if one identifies momentum with mu l (see below). 

To summarize our procedure, we suppose that equations (11) are correct (apart from 
higher multipole terms of course) and since ultimately, they were derived from (1), they 
have to be covariant. This claim leads to a constraint, the spin evolution equation. There 



/pZ pZ _i_ pZ pn pZ pn \ c k i 

2 ^ ki,m km,i ' nm ki ni kmJ^l u 

I 1 (-pi pn -pi -pn \ q hi _i_ ^ pZ ok 

O nm ki L ni L km) D l U ~T~ 9 km °l ■ 



. Thus, we 
(11) 
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is a certain similarity between this method and the original derivation of Papapetrou j2] . 
Papapetrou begins with a symmetric stress-energy tensor and finds, after having expanded 
up to the dipole order, an expression for the integrated stress-energy tensor that is not 
apparently symmetric. The claim that it should be symmetric then leads to a constraint 
which is just the spin evolution equation. 

Last, we remark that there is a tiny difference between (10) and (12) and the equations 
originally derived by Papapetrou in the sense that the momentum P l of the latter is 
identified with mu % in our equations, whereas the correct relation reads Pi = mui + 
DSikU k . We will derive the correct equations in the next section; we note however that 
the difference can be shown to be of second order in the spin and of first order in the 
curvature (see |3] or [T2|) and is thus negligibly small compared to the other terms. On 
the other hand, it has been shown in ^3] that (10) and (12) can actually be derived from 
the original Papapetrou equations by a slight redefinition of the center of mass of the 
body and can thus equally well be considered as exact (since we have not specified our 
definition of center of mass). (See also [5] on the issue of the center of mass redefinition.) 

3 The pole-dipole particle in Riemann-Cartan space 

We now turn to a classical pole-dipole particle moving in Riemann-Cartan geometry. We 
know that, since torsion couples only to the intrinsic spin, the classical pole-dipole particle 
should behave exactly in the same way as in a purely Riemannian space. It is however 
necessary to have at our disposal a consistent method of re-deriving the same equations in 
the full Riemann-Cartan framework for the later generalization to particles with intrinsic 
spin. 

Let us first give a short review of the basic concepts of Riemann-Cartan geometry 
and fix our notations and conventions. For a complete introduction into the subject, the 
reader may consult ^1]-P7] and the extensive reference list in [T7] . 

Latin letters from the beginning of the alphabet (a, b,c. . .) run from to 3 and are 
(flat) tangent space indices. Especially, rj a b is the Minkowski metric diag(l, — 1, — 1, — 1) 
in tangent space. Latin letters from the middle of the alphabet (i,j, k . . .) are indices in 
a curved spacetime with metric as before. We introduce the independent gauge fields, 
the tetrad e a m and the connection T ab m (antisymmetric in ab), as well as the corresponding 
field strengths, the curvature and torsion tensors 

ab -pab -pab _i_ -pa -pcb -pa -pcb 

Im m,l l,m ' cl m cm I 

rpa a _ a , b pa _ p opa 

Im m,l l,m ~ m x bl W- 1 - bm- 
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The spacetime connection Y % ml and the spacetime metric can now be defined through 



a _,_ pa b ap 



2 

mi 



It is understood that there exists an inverse to the tetrad, such that efe l b — 5 b . It can 
easily be shown that the connection splits into two parts, 



-pab pa6 , jy'db 

m m ' mi 



such that f ab m is torsion-free and is essentially a function of e^. K ab m is the contortion 
tensor (see below). Especially, the spacetime connection T l ml constructed from 



a _|_ pa b _ ap 



mi 



is just the (symmetric) Christoffel connection of general relativity, a function of the metric 
only. This notation is consistent with that of the previous section. 

The gauge fields e a m and T ab m are covector fields with respect to the spacetime index 
m. Under a local Lorentz transformation in tangent space, A a b (x m ), they transform as 



a A a b pa > A a A dye _ a a a c ho\ 



The torsion and curvature are Lorentz tensors with respect to their tangent space indices 
as is easily shown. The contortion K ab m is also a Lorentz tensor and is related to the 
torsion through K\ m = l m + T m % - T i , J, with K\ m = e i a e lh K ab m and analogously 
for T\ m . The inverse relation is T\ m = —2K l ^ m y 

Although we introduce the torsion tensor, which can be interpreted as the field strength 
corresponding to the translational group, we do not consider (local) translations here. We 
require our theory to be invariant with respect to local (tangent space) Lorentz transfor- 
mations and general coordinate (spacetime) transformations only. The identification of 
the tetrad (or its nontrivial part) as translational gauge potential needs a special treat- 
ment and is not needed here. (See however section 8.) 

All quantities constructed from the torsion-free connection T ab m or Y\ m will be denoted 
with a hat, for instance R d km = e t a e l b R ab km is the usual Riemann curvature tensor. Fur- 
thermore, we will use covariant derivatives with respect to proper time along a trajectory 
in the form 

c\P 

Dp - arm _ p p. v k 

ur ra — ^ 1 mk r i a 

DP a = ^-T b ak P b u k . 
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It is easily shown that DP m = e^DF s , which justifies the use of the same symbol D. Of 
course, the same derivative when referring to the torsion-less connection will be denoted 
byD. 

For the sake of completeness, and to avoid confusion concerning sign conventions, let 
us also derive the conservation laws for the stress-energy tensor and the intrinsic spin 
density (for details, see [14j). The stress-energy and spin density are defined by variation 
of the matter Lagrangian density C m as follows 

T ° = h (14) 

= esr»- (15) 

Under an infinitesimal Lorentz transformation (13), with A a b = 5% + e a b , where e ab = —e ba , 
the fields transform as 

SY ab = -e ab - T a e cb - T b e ac 5e m = e c e m (16) 

Therefore, the variation of the matter Lagrangian reads, omitting a total divergence 

6C m = 6 -^6e™ + ^ST ab m = e(2TH + D m a acm )e m (17) 

where D m is defined to act with T ab m on tangent space indices, and with T l ki (torsion free) 
on spacetime indices. The requirement of Lorentz invariance, i.e. SC m = 0, then leads to 

2T [ac] + B m a acm = 0. (18) 

On the other hand, the connection T ab m and the tetrad field e a m are covariant spacetime 
vectors with respect to m (thus, is a contravariant vector), and under an infinitesimal 
coordinate transformation 

x* = + (19) 

the fields transform as 

C(50=C(*) + ^, r ab m (x) = r ab m (x)-e, m r ab k . (20) 

We are interested in the change of the Lagrangian under active transformations. In order 
to evaluate the change of the fields at the same point x, we have to express the transformed 
fields with the old coordinates, i.e. 

= C(^) - O*)^ r a U*) = f ab m (x) - er ab m Ax). (21) 
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We then find for the variation, to first order in £ 

5e™ = ^ - (22) 

5T ab m = T ab m (x)-T ab m (x) = -e, m T ab k -er ab m , k . (23) 

The change in the Lagrangian therefore reads 

uc -a UL m 

= -2(eT a m e k a ), k C ~ 2eT a m e™ k e + (ea ab m T at \) ^ ~ ea ab m T ab m ^ k . (24) 
Requiring 5L m = and regrouping carefully the terms, we finally get 

(D m T m b )e b k + T m b T b mk = \R ab mk a ab m + \v ab k (V m a ab m + 2T [ab] ). (25) 

The last term vanishes with (18). The first term, written in terms of the usual (torsionless) 
covariant derivatives, (denoted with a semicolon ; ) reads 



r\ rpm b _ rpm b _ TZ-a rpm b 
U m-L b e k — 1 b;m e k 21 brn 1 a e k 

rpm jy-l rpm 

— 1 k;m 21 km 1 I J 

and the second term, using the relation T\ m = —2K l ^ lm ^ 

rpm rpb r r m (Tf^ T<T^ \ 

1 b 1 mk — 1 /I 21 mk 21 km)- 

Therefore, the first two terms of (25) simplify to T™. m — T™K l ' mk and the equation of 
motion finally takes the form 

1 
2 



rpm Ifi'mrp p>ab m (0P.\ 

1 k;m 21 fe 2 mi — mk° ah ■ \ ZD ) 



Recall that K tm k is antisymmetric in im, so that classical matter (with Tj mi ] = and 
a ab m = 0) will follow the general relativity relation T km m = 0. 

Let us now return to the classical pole dipole particle without intrinsic spin. We will 
re-derive the equations of the previous section using the concepts of Riemann-Cartan 
geometry. 

Based on the form of equation (8), we will try the following Lagrangian 

L = e\P a u { - h ab m S ah u m . (27) 
10 



Of course, we have coupled the (classical) spin to the torsion-less connection only. In (27), 
the momentum vector P a and the spin tensor S ab are considered as parameters and are 
neither functions of the coordinates nor of the velocities. They can be considered as the 
charges of the corresponding gauge fields e a m and T ab m - 

The Euler-Lagrange equations are then derived in a straightforward manner: 

p a „.i p p o. ip _ (-pab pab \ q i a p _ }^pab q (9%\ 

c i,m a r a c m,i a r a 2^ *i m m,i)°ab u ' e m r a ^ L m °ab- \^°) 

We use the same trick as before and express P a with the help of DP a and S a b with I)S a b. 
The result is 

DP m - T a m yP a + K a im u l P a = -±K* mi S ab u i + ab JI)S ab - u b P a + u a P b }. (29) 

It is easily shown, using DP m = DP m — K l mi Piu l and T l mi = -2^,,, that the left hand 
side is just DP m . 

In order for the equation to be Lorentz gauge invariant, we have to require 

DS ik = u k Pi - UiP k , (30) 
and then the momentum equation reduces to 

DP m = - l -R l \ m S lk u m . (31) 

If we define the mass through PiU 1 = m and require UiU 1 = 1, we can derive from (30) the 
relation 

P = mui + T)S ik u k . (32) 

Equations (30) and (31), together with (32) are exactly the equations derived by Papa- 
petrou in j2]. The equations derived in the previous paragraph differ by the second term 
in (32), which, as we have said, is of order S 2 . 

The derivation used in this paragraph shows clearly that the same equations are valid 
in every type of Riemann-Cartan spacetime. In the final equations, there is no torsion 
involved. The same equations hold, for instance, in a teleparallel theory. Torsion effects 
will only arise if we consider particles with intrinsic spin. 

4 Particles with intrinsic spin 

Particles with intrinsic spin require a more careful treatment, because one can easily run 
into problems, as has been clarified in In this section, we have in mind elementary 
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particles (extended bodies with intrinsic spin will be treated in the next section). There- 
fore, it seems, at first sight, obvious, to treat them as monopole particles, without dipole 
and higher order moments. Let us briefly review the Papapetrou method, as it has been 
applied in a Riemann-Cartan framework in [3] and [3], neglecting higher order poles. 

We start from equations (18) and (26). It is convenient to write them with spacetime 
indices only, and to use the following form 

(ea ) >m = -lei 1 J - 1 lm ea - 1 lm ea (33) 
( e ^"fc),jn — r^, m eT7 — K m k eT mi = -^R a mk <J u m , (34) 

where e = det = y/—g. Then, one considers a worldline X m and develops all the fields 
around this worldline, x m = X m + 5x m . Apart from (33) and (34), consider the following 
equations 

{ea lkm x n ), m = ea lkn + x n (ea tkm ), m (35) 
(eT>% = eT\ + x\eT m k ), m . (36) 

Integrating over three dimensional space and neglecting integrals containing Sx m (monopole 
approximation), one derives the following equations from (35) and (36) (for details, see 

0): 

u i p k = u ° J eT\d 3 x, with P k = J eT° k d 3 x, (37) 

a lm u k =u f e(T lrnk d 3 x ^ ^ ^k = f ^jiM^ (gg) 



where u % = dX*/dr. It is easy to show that a lk and Pk as defined above are tensors in 
the monopole approximation. Using those relations in (33) and (34), we arrive at the 
following equations 

Ba ik = (P i u k -P k u i ) (39) 
DPk-K^UrnPi = ^R u mk a u u m . (40) 

This completes the Papapetrou analysis of the monopole particle with intrinsic spin. 
However, as has been clarified in |4J, there are problems with the above approach. Indeed, 
the relation (38) cannot be considered to be generally valid. This is easily seen by con- 
sidering the Dirac particle, where a lkl is totally antisymmetric. For such a spin density, 
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the relation (38) leads to <j lk = 0, as can be seen by setting one index to zero and using 
the antisymmetry properties of a tkl . 

Applying WKB methods on Lagrangians for particles with integer and half-integer 
spins, the following relation that replaces (38) was derived in ^H] 



for a particle with spin s. This reduces to (38) in the limit of large spin, which coincides 
with the usual result for the Weyssenhoff spin fluid. Thus, (38) can be considered to hold 
for macroscopic spin polarized bodies, whereas elementary particles obey (41). For spin 
one-half, we get from (41) a totally antisymmetric spin density. 

In order to avoid the relation (38), that arises in the monopole approximation of 
the Papapetrou approach, the authors of |£] concluded that one has to include dipole 
moments of the particle, i.e. to introduce a classical spin besides the intrinsic spin. 
The same has been done in [5J, where the Mathisson-Papapetrou equations in Riemann- 
Cartan space were derived for the first time. This seems a rather strange concept for an 
elementary particle and is problematic from a practical point of view. How can one discern 
experimentally the spin and the rotational momentum of an elementary particle? How can 
we define the intrinsic spin of an electron if it is not the total angular momentum in the 
particle's rest frame? Therefore, throughout this article, we take the point of view that 
the intrinsic spin of an elementary particle corresponds to its total angular momentum 
in the rest frame of the particle. There are possible objections to this, especially when 
one considers compound particles, like the proton for instance. In such cases, one could 
say that the spin is partly due to the spin of the constituent quarks and partly of orbital 
nature (some kind of rotation of the quarks around the center of mass of the proton, in a 
semiclassical picture) . Indeed, experiments are carried out that aim at determining these 
different parts. However, the fact that the total spin appears always to be exactly 1/2 
leads us to believe that the proton as a whole can be treated quantum mechanically as 
a spin 1/2 particle. If the spin is partly of non- intrinsic nature, one has to answer the 
question why the orbital part cannot be transferred to other particles and why it always 
sums up with the intrinsic part to exactly 1/2. The final answer to this question, however, 
should belong to the experiment. To this aim, the full equations, containing intrinsic spin 
as well as dipole correction terms, as derived in jl] and [3], could eventually be used. 

In order to avoid the problem with equation (38) without taking into account dipole 
terms, we will choose another way of dealing with (semiclassical) elementary particles. 
Having in mind that the connection T ab m = T ab m +K ab m couples linearly to the spin density 
in elementary particle Lagrangians, we conclude from the fact that the spin density is of 




(41) 
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the form (41), that the spin-torsion coupling will be of the form 

K^laikU 1 + —a l iU k + —a k l Ui] = K* ik l( T ik u l . (42) 
2s 2s 

The left hand side determines K* th l to which we will refer as the effective torsion 

K *ik i = K ik i + l_ Ki ik + l_ K ki (43) 

Only this part of the torsion will couple to the spin. For a macroscopic spin-polarized 
body, K* tk l reduces to the full K lk l and for the Dirac particle (s = 1/2), it reduces to the 
totally antisymmetric part of the torsion. All quantities formed from the effective torsion 
will be denoted with a star, especially T* a ^ = T ab m + K* ab m . 

Let us now derive the equations of motion for spin and position for a point particle 
with intrinsic spin. If we suppose that in a purely Riemannian space, the particle with 
intrinsic spin will behave just like a particle with classical spin, we have to introduce 
the term — \T ab m a ab u m into our Lagrangian (see (27)). If we add to this the spin-torsion 
coupling (42), we are led to the following Lagrangian 

L = elPy - \T™ b m a ab u m . (44) 

The ultimate physical justification for the use of a specific Lagrangian can only be the 
correctness of the equations of motion derived from it. 

Just as in the previous sections, we readily derive the Euler-Lagrange equations in the 
form 

D*P m - T* a mt u l P a = -±R* ab mi <j ab v; + ^r* ab m [D*a ab - P a u b + P b u a ). (45) 
Requiring the gauge invariance under a local Lorentz transformation, we finally get 

B*P m -TZXPa = -\w ab mi a ab u\ (46) 
DV ife = PiU k - P k Ui. (47) 

In order to compare with the classical geodesies, we can rewrite the first equation as 

DP m - ijtSJPV - PV) = - l -R*^a ab u\ (48) 

For a macroscopic spin-polarized body, these equations agree completely with the result 
from the monopole approximation in j3], i.e. with equations (39) and (40). However, 
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such a macroscopic body will also possess dipole moments and the equations have to be 
adopted to that case (see next section). For the spin half particle (totally antisymmetric 
torsion), the results agree, if we ignore again the quantity P^u k \ which is of order a 2 , 
with the results from the WKB analysis of the Dirac equation, as carried out in [T^]. Also, 
in the same approximation, there is an agreement with the results of [2U| for the spin 3/2 
particle and of j2I] and (221 f° r the spin 1 or Proca particle. 

Thus, there is no need to consider dipole moments while dealing with elementary 
particles. The Lagrangian (44) provides a very comfortable method of deriving the correct 
equations of motion without reference neither to the gravitational field equations, nor to 
the specific matter Lagrangian for the particle under consideration. The only additional 
information we needed was relation (41). 

Equations (46) and (47) were derived in |7j for the case of a purely Riemannian space 
(T a lm = 0), following a slightly different procedure. The term T ab m a a bU m was interpreted 
as the Hamiltonian part of a Routhian R and the spin evolution equations were then 
derived using the Heisenberg equations = i[R, er &] and the commutation relations 
for the spin tensor. However, in order to get the correct equations, they had to redefine 
the spin tensor and to modify slightly the Routhian. This was due to the fact that their 
Lagrangian part was taken to be of the form myfuiU 1 . On the other hand, if we start 
with the form e°;P a u l and suppose that P a and <r a b obey the commutation relations of the 
Poincare algebra, one gets the correct precession equation following the same procedure 
as in [7j. We can see this as an independent confirmation of our method, since the mere 
requirement of Lorentz gauge invariance leads to the same result. 

In 0, a different method was used, involving additional variables to describe the spin 
degrees of freedom. However, although initially intended to describe elementary particles 
in Riemann-Cartan spacetimes, the equations actually correspond to the special case 
s — ► oo. The same holds true for the equations given in [2B| and [21]. Another, more 
recent attempt to generalize the multipole formalism to Riemann-Cartan geometry can 
be found in [25] . 

5 Macroscopic bodies with intrinsic spin 

Finally, we are ready to consider macroscopic spin-polarized bodies. This is probably the 
experimentally most important case, since gravitational effects on elementary particles are 
usually very small. (See however [2H|-[2E] for experiments involving gravity on a quantum 
mechanical scale.) In this section, we have in mind an extended body, like a neutron 
star, for which the spins of the constituent particles are (partially or fully) aligned. In 
addition, the body may rotate. From the last two sections, we may expect a coupling of 
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the classical spin S ab to the Riemannian connection f ab m , and a coupling of the intrinsic 
spin a ab to the effective connection, which coincides in the macroscopic case with the full 
connection r* a m = T ab m . Hence we start with 

L = elPaU* - ~f ab m S ab u m - l -T ab m a ab u m . (49) 

We derive the following equations 

I)P m - l -K lim {P l u l - P l u l ) = -^R lk mi S lk u* - l -R ah mi a ab u l 

+ l -T ab jDa ab - P a u b + P b u a ] + if ab m BS ab . (50) 

We can now regroup all the terms proportional to T ab m , or equivalently to T ab m , and then 
require Lorentz covariance. The result is 

DP m - X -K lim {P l n l - P l u l ) = - l -R lk mi S lk u* - ^R ab mi a ab u l - ^K* k m DS ik (51) 
and for the spin evolution 

Da ik + DS ik = PiU k - P k Ui. (52) 

We can use (52) in (51) to eliminate the apparent coupling of the torsion to the rotational 
moment in the last term. The equations are more or less what could have been expected 
right from the start, except maybe for the last term in (51). Our equations agree with 
those derived in jl] (for the macroscopic limit s — > oo) using the Papapetrou method in 
the dipole approximation. 

These equations alone are of course not sufficient to determine the behavior of the body. 
An additional relation between the internal and the classical spin has to be assumed. For 
a simple model of a neutron star, let us suppose a strong coupling of the form Si k = ga ik 
with constant g. In this case, we can write (52) in the form 
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PiU k - P k Ui 



Now, if we introduce the connection 



T/m — Tim + ^ _|_ ^\m) (53) 



we can write ^ 

t)S ik = PiU k - P k Ui. 
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Multiplying with u 1 , denning the mass as m = PiU 1 and using 1 = UiU 1 , we find 

P k = mu k + -— - DS ki u\ (54) 
9 

and the precession equation finally reads 

I)S ik = (u k bS u - UiDS kl )u l . (55) 

Thus, the intrinsic spin (and also the classical as well as the total spin) is Fermi- Walker 
transported with respect to the connection (53). 

Introducing the constraint Si k = gOi k directly into (49), we see that the connection 
(32) governs the whole propagation, and not only the spin evolution. The equation for 
momentum propagation is easily derived then. 

We used the constraint Si k = g<Jik mainly as an illustrative example. For realistic 
neutron star models, you can consult (2H1 and (301 for instance. You should however have 
in mind that these models rely strongly on general relativity as underlying gravitational 
theory. Since the use of our equations of motion naturally suppose that the underlying 
theory is a Poincare gauge theory with dynamical torsion fields, these models have to 
be revisited in view of the modified gravitational interaction. This might lead to severe 
changes, since the interior of neutron stars is governed by very strong fields. To derive 
such models is beyond the scope of this article. 



6 The Dirac particle 

In order to gain more confidence in the procedure we followed in section 4, mainly the in 
introduction of the effective torsion (43), we will have a closer look at the most important 
example, that of a spin 1/2 particle and present an alternative derivation of the equations 
of motion. 

The Dirac particle is described by the following Lagrangian density 

C= l - [^ m D m ^ - D m ^ m ^] + m^, (56) 

where D m ^ = (d m - ^T ab m a ab )^j, D m ^j = + ^T^a ab and 7 m = e™ 7 a , with 7 a the 
usual Dirac matrices and o ah = ^[j a , 7b]- The Dirac equation reads (T m = T* mi ) 

i(j m D m -j m T m )Tp = rmp. (57) 
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The stress-energy tensor (2e) 1 5(eC)/8e l a is found to be 

^ = 1$7°I^ - A^7°# (58) 

Several remarks are in order at this point. It is a well known fact that only the totally 
antisymmetric part of the torsion couples to the Dirac particle. This follows immediately 
from (57), which can be written as 

i! m D* m i, = (59) 

where D* m is constructed from the connection (43) with s = 1/2. Furthermore, the 
Lagrangian (56) is numerically equal to 

C* = % - [^ m D* m ^ - P>^ 7 "V] + m^. (60) 

This Lagrangian is preferable from the point of view that it incorporates a minimal cou- 
pling principle that can be applied equally well in the Lagrangian or in the field equations. 
However, (56) and (60), although they both lead to the same Dirac equation, are not com- 
pletely equivalent, since for the latter the stress-energy tensor becomes 

r ! = ^rt-w# (6i) 

which differs from (58) in the non-axial torsion parts contained in the latter. Once again 
(61) seems preferable, since the Dirac particle does not couple to these torsion parts. 
However, it makes of course a difference which of these tensors actually represents the 
source term of the gravitational equations. If we remain in the logic of a Lorentz gauge 
theory and have in mind that there might be additional sources with different spin (cou- 
pling to other torsion parts), the only consistent way is to write down only Lagrangians 
that contain the full connection T ab m , as is the case for (56), even though they might 
contain parts which do not couple to the particle in question. Otherwise we can not carry 
out the variation with respect to T ab m . 

However, as long as we are interested only in the equations of motion for the Dirac 
particle (or their semiclassical limit), and not in the gravitational field equations, the 
use of either T ab m or T* ab m is equivalent. Equations (57) and (59) are exactly the same. 
Under this aspect, the use of T* ab m is certainly preferable. It can be seen as the physically 
relevant connection since it is the connection that really couples to the Dirac particle. 

What does this have to do with our equations of motion? Well, our point is that, if we 
start, as in the usual Papapetrou procedure, from the divergence relations for the stress- 
energy tensor and the spin density, and then define the momentum as Pj = f eT^d 3 ^, 
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this is already very unfortunate right from the start, since this momentum vector contains 
implicitly non-axial torsion parts from (58) which should not couple to the Dirac particle 
at all. This explains at least partly why the monopole approximation does not lead to 
correct results for elementary particles A first step should be to look for a divergence 
relation involving the tensor T*" and carry out the Papapetrou method on this tensor. 

The above considerations are also useful in another aspect. The use of the covariant 
derivative D* m allows us to generalize the proper time formalism, i.e. the covariant gen- 
eralization of the Heisenberg equations, to the case of a Riemann-Cartan space. We will 
briefly discuss the concepts and refer to and the references therein for further 

details. Notably, in .34], the resulting equations were used to extract for the first time 
the precession frequency of the spin of a Dirac particle in an axial torsion field. 

One starts with the operator H whose eigenvalue is (the negative of) the mass, Hip = 
—mip. (The mass is the only available covariant generalization of the classical energy, 
apart from the stress-energy tensor.) Then identify (for any operator A) 

A A 

i[H,A} = — , (62) 

where r is interpreted as proper time along the semiclassical trajectory. In flat spacetime, 
the Hamiltonian is just H = —j l Pi, with pi = idi and EM fields can easily be taken into 
account by pi — > 7Tj = idi + eAj. With this in mind, and having a look at (59), we find 
that the Hamiltonian for our case has the form 



H = -i m P m , (63) 

with 

P m = iD* m = i{d m - l -Y ah *a ah ), (64) 

where T a6 ^ is again the connection which contains only the totally antisymmetric part of 
the torsion. For the velocity operator, we find 

u l = — = i[H,x 1 } =7*. (65) 
dr 

We thus have the expected relation H = —u m P m = —m. Let us have a look at the 
Heisenberg algebra for momentum and position operators. We have: 

[x\x k ] = 0, [P t) x k ]=i5t 
[Pi, P k ) = -R* c ik a cd . 
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This is in complete analogy with the EM case (just recall the relation [7^, ir k ] = ieF^). A 
useful relation is also [Pi,j k ] = -i^ l T* k H = -i^ l (T* k a + T* k u ). 

Writing down the commutation relations for P m ,7* and a a b with the Hamiltonian, it 
is straightforward to derive the following equations of motion 

DVjfe = P\u k - P k ui (66) 
V*P m = ~R* cd mk u k a cd + T* k ml u l P k (67) 
DV = -a mi P m . (68) 

Of course, D* appearing here is again the proper time derivative acting on tensors (not 
on spinors) with T*\ m . Relation (68) has been related to the zitterbewegung and has no 
classical counterpart (see [HI]). The other two are exactly our equations (46)-(47). 

Both derivations, the one of section 4 and the present one, although concept ionally 
quite different, lead to the correct results of the WKB analysis |19j . In both cases, the 
crucial step was to use the correct connection containing only the part of the torsion that 
effectively couples to the spin density. 



7 Spin-polarized test body in Brans-Dicke theory 

Finally, we will take a look at a specific example involving the equations derived in section 
5. Since in most practical cases, intrinsic spin effects will be very small (see however |26| - 
[2*S]). we turn to a theory that gives rise to a torsion field even in the absence of a spinning 
source, namely Brans-Dicke theory generalized to Riemann-Cartan geometry. The theory 
is based on the action 

S = J e Lr - ^±1 _ c \ d 4 x _ (6Q) 

It is easily shown by solving the equation that arises from variation with respect to the 
connection, that in the case of a spinless matter Lagrangian (the source), we get a torsion 
field of the form (cf. [SSI, OT> 

The field and the metric g^. are subject to the classical Brans-Dicke field equations, as 
can be seen by plugging (49) into the equations for <p and e^. The torsion (70) is entirely 
determined by its trace alone. It is often claimed (probably having the Dirac particle in 
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mind), that such a torsion field is pure gauge, and does not lead to physical consequences. 
However, if the spin of the test particle moving in this field is not totally antisymmetric, 
there might well be a measurable effect. 

We are interested in a macroscopic test body with intrinsic spin moving in the field 
of a classical, spherically symmetric source. (For equations of motion of a classical point 
mass in Brans-Dicke theory, see for instance |37| and references therein.) For simplicity, 
although not quite realistic, let us suppose that the test body does not rotate. The 
spin evolution is then described by equation (52) with S^. = 0. We proceed as in jlj 
and introduce the spin vector a 1 = jrf klm Ukcri m where rf klm = ls lklm . If we impose the 
condition <JkiU k = we can invert to a = —r] %klm Uia m . Omitting the higher order terms, 
we are left with 

Da* = DV + K\ m a l u m = 0. (71) 

We are not interested in the momentum equation, because the deviations from the 
geodesies are expected to be very small. 

Since the metric and scalar field are determined from the same equations as in classical 
Brans-Dicke theory, we can directly use the results from the post-Newtonian expansion 
as given in [3Bi for instance, and write 

H> 
9oo 

9a(3 

The coefficients read 

a 
b 

c 

Here, m is the mass of the source body and s its sensitivity (see [SB] and references 
therein). The sensitivity is basically the binding energy per unit mass. Recall that 
the static spherically symmetric solution of Brans-Dicke theory is not a one parameter 
solution as in general relativity, but depends on more properties of the central source. In 
the parameterized post-Newtonian formalism, it turned out that an astrophysical body 
can conveniently be described by its mass and its sensitivity. The latter takes values from 
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10 6 for a usual star (like the sun) to one half for black holes. Neutron stars have values 
around 0.1 - 0.3 (see [SHI). 

We now introduce (72)-(74) into (70) and (71). From the construction of cr\ we have 
a l Ui = 0. Using this relation and parameterizing the equation with the time coordinate, 
we get the result 



da m 
dt r 3 



a + 26 + c,=> _ 
(L x cr) 



a c 
— 2^(^ ' + ' %)) + ba(x ■ v) + -(v(x ■ a) + x (v ■ a)) , (78) 

with L = v x x the orbital angular momentum per unit mass of the test body. We have 
separated terms antisymmetric and symmetric in (x, v) . Just as in general relativity, the 
latter can easily be shown to have mean values over one orbit that vanish for closed orbits 
and that are negligible for the quasi-closed orbits we consider here (i.e. orbits closed, up 
to a small spin precession correction). 

Thus, the relevant precession equation reads: 

d " m a + 2b > —Lxa). (79) 



dt r 3 2 

The torsion corrections come from the term in c. 

In the limit u — > oo, the parameters take the values a = l,6=l,c = and we get the 
well known factor | from general relativity in the spin precession equation, but this time 
for intrinsic spin. This confirms once again the result that in general relativity, a particle 
with intrinsic spin behaves just like a particle with classical spin. More precisely, in the 
u — > oo limit, the theory goes over not to general relativity, but to Einstein-Cartan theory 
(see jHHj for a recent review). Since we consider only spinless sources however, there is no 
difference in the resulting solutions. (Even for a spinning source, the differences vanish 
outside the source.) 

In Brans-Dicke theory, the only case where the torsion vanishes is when s = |, i.e. 
when the source is a black hole. In this case, c = and the scalar field is constant. 

The interesting cases are of course when c ^ 0. The largest influence of the torsion 
is found for ordinary stars, when s ~ 0. We can evaluate for this case the ratio of the 
torsional part and the other parts for different values of u, for instance 

cj^ 500 -> ^--^1,4-10~ 3 



a + 26 

c 

a + 2b 



uj w 1000 -> w 6 • 10~ 4 . 
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These results can be interpreted in two ways: Within the framework of Brans-Dicke theory 
with torsion, they tell us the difference in the precession equation of a classical rotating 
test body (which does not couple to torsion) and a spin polarized test body. Thus, we 
can in principle determine, whether a body is just rotating or possesses intrinsic spin. 

On the other hand, for a spin-polarized body, they tell us the difference between the 
precession in classical Brans-Dicke theory without torsion for which the spin precession is 
the same for rotational and intrinsic spin and the precession of the same body in Brans- 
Dicke theory in a Riemann-Cartan geometry. 

The results are easily generalized to the more realistic case of a rotating test body 
using (55). More work is involved to get results for the two body problem (for instance 
a binary consisting of two rotating, spin polarized neutron stars). However, as far as the 
test body is concerned, we have shown that already the pure spin terms are very small. 
If there is an additional rotation, the relevant connection will be of the form (53), where 
for a realistic neutron star model (Sik » crik), we wm have a large g and therefore, the 
torsion effects will become even smaller. 

As far as the source is concerned, in a realistic neutron star model, the rotational 
effects will by far dominate the spin effects, so that the field will essentially be a Kerr-like 
analogue of Brans-Dicke theory. In any case, the torsion outside the source will again be 
produced by the scalar field only, because the spin of the source gives rise only to non- 
dynamical torsion fields (vanishing outside of the spin density), just as in Einstein-Cartan 
theory. A discussion of geodesies and autoparallels (for classical point particles) in such a 
Kerr-Brans-Dicke geometry can be found in 0U|. These results may easily be generalized 
to take into account the spin effects. 

Nevertheless, the analysis shows that there are physical consequences of the vector 
torsion, although very small, and we cannot consider it as a pure gauge field. 

8 Symmetries of the Lagrangian 

Throughout this paper, we have obtained the spin evolution equation by requiring the 
momentum equation that results from the Euler-Lagrange equations of the Lagrangian 
under investigation to be either covariant under spacetime diffeomorphisms in the Rie- 
mannian case or under local Lorentz transformations in the Riemann-Cartan case. One 
might wonder whether the same argument can be applied directly to the Lagrangian, i.e. 
if the requirement of the invariance of the Lagrangian under those transformations leads 
to the equation for the spin evolution. We will investigate this subject case by case. 
Let us begin with the Lagrangian (8). The only symmetry that occurs in Riemannian 
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geometry is the diffeomorphism invariance, which can be written in the infinitesimal form 

x i ->z i + e < (z). (80) 

This contains Lorentz transformations (e l (x) = EikX k with antisymmetric and con- 
stant) as well as translations (e l = a 1 with constant a 1 ) as subcases. The requirement of 
the Lagrangian to be invariant (up to a proper time derivative) under (80) leads to the 
Euler-Lagrange equations if the Lagrangian is a function of x l and u l alone. However, in 
the case of (8), there are additional quantities 5^ that transform as a tensor under coordi- 
nate transformations. Therefore, the requirement 5L = leads to the equation 5^ = 0, as 
is easily shown. This cannot be used as a constraint; it is neither covariant, nor physically 
acceptable as equation of motion. Therefore, we can conclude that the Lagrangian (8) is 
not a scalar. Covariant under (80) are only the final equations of motions. 

However, in the case of the Lagrangians (27), (44) and (49), no such problems occur, 
since the quantities P a and S a b transform as scalars under coordinate transformations. 
The requirement for L to be invariant under (80) simply leads to the Euler-Lagrange 
equations. In those cases, the scalar character of the Lagrangians is obvious anyway, 
since the connection as well as the tetrad field transform as covariant spacetime vectors. 

Let us turn to the gauge transformations. We consider the Lagrangian 

L = e^Pj - h ab m S ab u m , (81) 

which stands exemplary for the Lagrangians (27), (44) or (49). (Remember that the 
connections, T ab ml Y ab m and T* a ^ all transform in the same way.) Before turning to the 
Lorentz transformations, let us make a remark on gauge translations. 

In Poincare gauge theory, the ten gauge fields T ab m (Lorentz) and T a m (translations) 
alone are not sufficient to write down an invariant Lagrangian. This is due to the fact that 
the translational field T a m does not transform homogeneously as Lorentz vector. (It might 
even vanish for instance and is certainly not invertible, and thus not suitable to be used 
as a tetrad field.) One therefore introduces an additional field £ a (the coset parameters) 
transforming as £ a — > £ a + e a b !; b + e a under Poincare transformations and defines the tetrad 
field as e a m = T a m + D m e a . Then, this tetrad field transforms homogeneously as Lorentz 
vector (i.e. it is invariant under translations). Therefore, every Lagrangian expressed in 
terms of the fields T ab m and is trivially invariant under translations. The translational 
gauge symmetry is hidden. From this point, we started our discussion of Riemann-Cartan 
geometry in section 3. The details and the foundation of this interpretation of Poincare 
gauge theory can be found in [4*Tj . 

This leaves us with the local Lorentz transformations. The field transformations as 
given in section 3 can be written for infinitesimal Lorentz transformations with parameters 
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e ab (x) = —e ba in the form 

5T ab m = -D m s ab and 5e a m = e\e h m . (82) 

Let us for one moment assume that the quantities P a and S ab in (81) are invariant (i.e. 
scalars). Then, the change in the Lagrangian reads 

5L = e\e b m P a u m + l -Y) m e ab S ab . (83) 

Dropping a total proper time derivative and using the antisymmetry of e ab , this can be 
written in the form 

SL = -\e a \VS ab - (P a u b - P b u a )). (84) 
Thus, the requirement 5L = leads to the correct spin evolution equation DS^ = P a u b — 

PbU a - 

However, this result was derived assuming that S ab and P a are invariant under Lorentz 
transformations. Since they are ultimately related to spacetime tensors via = e1e\S ab 
and Pi = e^Pa, we know that they have to transform as Lorentz tensors, as has been 
assumed throughout the article. Including this additional variation in (83) leads to 5L = 
— \e ab S a b- If we require this to vanish, we are let to equation S ab = 0, which not only is 
not covariant, but has also to be rejected on physical grounds. 

Therefore, we conclude that the Lagrangians (27), (44) and (49) are not invariant 
under local Lorentz gauge transformations. Covariant are only the final equations of 
motion. 

9 Conclusion 

Our main result is summarized in the Lagrangians (44) for elementary particles and (49) 
for macroscopic bodies. Varying with respect to the coordinates and requiring gauge or 
diffeomorphism invariance of the resulting equation, the equations of motion for spin and 
position are obtained very easily and without any ambiguities. 

The method can be justified by the fact that the results are in agreement with the WKB 
analysis of elementary particles and with those obtained from the pole-dipole expansion 
following Papapetrou's method. 

We did not discuss in general the final equations in this paper, because they are not 
new, except for the result that a rotating, spin-polarized body with a strong spin-rotation 
coupling S = ga, couples to the connection given by (53). Nevertheless, we took the 
opportunity to show briefly that the vector-torsion arising in generalized Brans-Dicke 
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theory affects the spin precession of a macroscopic body and should not be referred to as 
pure gauge. 

Finally, we showed that the Lagrangians under investigation are not scalars under the 
transformations under consideration, nor does the requirement for them to be invariant 
lead to the correct spin evolution equations. This means that the processes of varying 
with respect to the coordinates and requiring gauge or diffeomorphism invariance do not 
commute. 
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